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- - - Here we examine various aspects of the pairing mechanism for two models, the two-dimensional and two- 

leg ladder Hubbard models. Both of these models exhibit pairing correlations with d 2_y2 symmetry. 
However, the undoped insulating states of these two systems are different with the two-dimensional lattice 
characterized by a ground state with long-range antiferromagnetic order and the two-leg ladder having a 
spin gap in its ground state. Our aim is to gain a better understanding of the underlying pairing mechanism 
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which causes d^2_y2 pairing in these two models. 



OO ■ 1 Introduction 

^ . 

Understanding the pairing mechanism responsible for superconductivity in the high-temperature 
^ , superconducting cuprates remains one of the central issues in condensed matter physics. Various 

OO ' experiments strongly suggest that the dominant symmetry of the gap in these materials is 

, dx2-y2. However, the implications of this for the pairing mechanism remain controversial |^]. In 

' one view , antiferromagnetic spin- fluctuation exchange models provide a simple framework 

, for understanding the dx2-y2 symmetry of the gap. In addition, various phenomcnological cal- 

\^ ' culations based upon these ideas have often provided remarkable fits to data However, the 

0> , effect of vertex corrections H as well as the unusual normal state properties and in particular 

the pseudogap phenomena ||9|] observed in the under-doped materials raise key questions which 
. , suggest that important pieces of the puzzle remain to be understood. Thus it seems useful to 

Vh ' review what is understood about systems which show d^2_y2 pairing fluctuations. Specifically, we 

propose to review some results which have been obtained for the two-dimensional Hubbard model 
^ ' and the two-leg Hubbard ladder. While it is not known whether the doped two-dimensional Hub- 

O . bard model has a low-temperature superconducting phase, there is numerical evidence for dx2-y2 

pairing fluctuations in this model mi^,|ll|. For the doped two-leg ladder, groundstate calcula- 
tions Jl2| clearly show the presence of (ia;2_j,2-like pairs. Furthermore, the undoped insulating 
states for these two systems are different. The ground state of the two-dimensional Hubbard 
model has long-range antiferromagnetic order jl^ , while that of the two-leg Hubbard ladder is 
■ characterized by a spin gap [Q. Both, however, have strong short-range antiferromagnetic cor- 

relations. Thus it is of interest to compare them and here we will examine these two cases with 
the aim of gaining insight into the mechanism underlying the tendency towards dx2_y2 pairing 
in these systems. 

After a summary of the energy scales which arise in the Hubbard model, we first turn to a 
discussion of the two-dimensional Hubbard model. Here we note the development of short-range 
antiferromagnetic correlations in the doped two-dimensional Hubbard model as the temperature 
is lowered below the exchange interaction At'^/U. We examine the effect of these correlations 
on the structure of the quasiparticle spectrum and the effective particle-particle interaction. We 
then discuss how this interaction leads to a tendency towards d^2_y2 pairing. 

Following this, we turn to the example of the two-leg Hubbard ladder. Here we first discuss 
the half-filled insulating case, distinguishing between a band insulator and a strongly correlated 
spin gap insulator. Then we study the case of a doped ladder and discuss the nature of the 
observed pairing correlations. 
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The Hubbard model provides a simple model of the Cu02 system: 

H ^-tY^ (4cjs + c]^Cis) n.in^. (1) 

Here cj^, creates an electron of spin s on site i, in the first term the sum is over near-neighbor sites, 
and Tiis = c|^Cis is the occupation number for electrons with spin s on site i. The one-electron 
transfer between near-neighbor sites is t, and U is an onsite Coulomb energy. In Eq. (1), we 
have chosen to keep only a near-neighbor hopping, but one could of course also add a next-near- 
neighbor hopping t' . 

The bare energy scale is set by the bandwidth 8 1 and the effective Coulomb interaction U, 
which are both of order several electron volts. Near half-filling, electrons on neighboring sites tend 
to align antiferromagnetically so as to lower their energy by the exchange interaction J = At^ /U . 
This interaction is of order a tenth of an electron volt and, as we will see, sets the energy, or 
temperature scale, below which antiferromagnetic correlations, the low-energy structure in the 
single-particle spectral weight, and the pairing interaction develop. 

While Monte Carlo and Lanczos calculations for a 4 x 4 lattice find that two holes 
added to the half-filled Hubbard ground state form a dx2_y2 bound state, and density matrix 
renormalization group calculations [ p2| find that dr^2 _y2-\\ke pairs are formed on two-leg Hubbard 
ladders, it is not known what happens for the two-dimensional Hubbard model. It is possible 
that on an energy scale of order J/10, a dx2_y2 superconducting state forms. However, this 
may well require modifications of the model, such as an additional near-neighbor A JSi • Sj term 
or possibly a next- near-neighbor hopping t' . Nevertheless, it is known that as the temperature 
is reduced below J, dj.2_y2 pairing correlations develop in the doped two-dimensional Hubbard 
model, and here we will examine why this happens. 

2 Two-Dimensional Hubbard Lattice 

At half-filling, {ni-^+riii) = 1, the 2D Hubbard model develops long-range antiferromagnetic order 
as the temperature goes to zero jl^. In the doped case, strong short-range antiferromagnetic 
correlations develop as the temperature decreases below J. This is clearly seen in the temperature 
dependence of the wave vector and Matsubara frequency-dependent magnetic susceptibility 

X(q,zc.™) ^j-Y, r dre^---e-^^-^{m-_^,{T)mt{0)) . (2) 

e ° 

Here m^{Q) — c|^Ci| and to^^(t) = e^'^m^_^g{0)e^^'^ , where rn~^g{0) is the hermitian conjugate 
of m^^(O). Monte Carlo results for x(QiO) versus q along the (1, 1) axis for an 8 x 8 lattice 
with U/t = A and a filling (n) — 0.875 are shown in Fig. |l|(a). In Fig. |l|(b), the Matsubara 
frequency dependence of versus tOm — 2nnrT is shown for q — (tt, tt). The inset of 

Fig. |l|(b) shows the temperature dependence of the antiferromagnetic correlation length Here, 

is defined as the half-width at half-maximum of x(ctjO). The low-frequency nature of the 
antiferromagnetic correlations are seen in Fig. ^(a) where Imx(q = (tt, tt), uj) versus io is plotted 
at various temperatures. In addition, Fig. ^(b) shows Imx(q, versus q in the ujz 
limit. This is the quantity which determines the NMR Tf^ response of the system. These figures 
clearly show the development of significant short-range and low-frequency antiferromagnetic 
correlations as the temperature decreases below J ^ At^/U. 

As these antiferromagnetic correlations develop, the single-particle spectral weight 

A{p, ui) — ——Im G(p, iujn LO + iS) (3) 
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Figure 1: (a) Magnetic susceptibility x(cii 0) versus q along the (1,1) direction for various temperatures. 
These results are for an 8 x 8 lattice with U/t = 4 and a filling (n) = 0.875. (b) Matsubara frequency 
dependence of x(ci, i^m) for q — (tt, tt) at the same temperatures as in Fig. 1(a). Note that as T decreases 
below 4t^ /U, short-range antiferromagnetic fluctuations develop. Inset: Temperature dependence of the 
antiferromagnetic correlation length ^. 



and the density of states 

p 

also change. Figure ||(a) shows the temperature evolution of N{uj) for U/t — 8 and (n) = 0.875. 
Figure||(b) shows N(u;) for U/t = 12 and T — 0.5t. As the temperature is lowered, a peak appears 
on the upper edge of the lower Hubbard band. This peak arises from a narrow quasiparticle band 
shown in the single-particle spectral weight A{p, lu) of Fig. ^ and plotted as the solid curve in 
Fig. ||. As the momentum p goes towards the F point (0,0), we believe that the quasiparticle 
peak is obscured by the lower Hubbard band because of the resolution of the maximum entropy 
technique which we have used. Indeed, at the F point a separate quasiparticle peak is found 
from Lanczos exact diagonalization on a 4 x 4 lattice. The solid line in Fig. || shows what 
we believe is the quasiparticle dispersion relation. As clearly evident in the spectral weight 
shown in Figs, ^(a) and (b), the quasiparticle dispersion is anomalously flat near the (tt, 0) 
corner. As discussed by various authors [[ITyiql , this reflects the influence of the antiferromagnetic 
correlations on the quasiparticle excitation energy. It is clear that the peak structure in N(uj) 
also arises from the short-range antiferromagnetic correlations and is a many-body effect rather 
than simply a non-interacting band Van Hove singularity. 

Monte Carlo calculations have also been used to determine the two-fermion scattering vertex 
r(pi,P2,P3,P4)- Here pi stands for (pi,ia;„J and a. Using F and the single-particle Green's 
function G{p), one can solve the i-matrix equations shown diagrammatically in Fig. |^ to obtain 
the irreducible particle-hole and particle-particle vertices. The singlet irreducible particle-particle 
vertex Tis{p' , —p',p, —p) in the zero center-of-mass momentum and energy channel represents the 
effective pairing interaction. In Fig. 0(a), Tis{q = p — p') is plotted for q along the (1, 1) direction 
and iujn — iuJn' = inT, corresponding to zero Matsubara energy transfer. The Matsubara 
frequency dependence of this vertex for energy transfer ujm = ujn' — ojn is shown in Fig. |^(b) 



3 



1 J 1 

; (a) 


1 1 1 1 1 1 
T=0.25t - 


T=0.33t ; 

---T=0.50t - 
--T=1.0t '. 


■ / V. 







3 

E 



5 
cj/t 



10 



5 - 




(0.0) (tt.O) , 







o. 



-Mr: 



(0.0) (tt.O) (tt.tt) 

q 



(0.0) 



Figure 2: (a) Spin-fluctuation spectral weight Im x(q, w) versus a; at q = (tt, tt) for various temperatures, 
(b) lmx{<i, <^z)/iOz versus q in the hmit of cjz — > for the same temperatures as in Fig. 2(a). Here q 
is plotted along the path shown in the inset. These figures are for U/t = 4 and (n) = 0.875. 




Figure 3: (a) Evolution of the single-particle density of states N{uj) with temperature for U/t = 8 and 
(n) = 0.875. (b) N{oj) versus oj for U/t = 12, (n) = 0.875 and T = 0.5t. The vertical dashed lines denote 
the chemical potential. 
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Figure 5: Dispersion of the quasiparticle peak in the spectral weight versus p for U/t = 8, (n) = 0.875 
and T — 0.5t. The sohd points mark the low-cncrgy peaks of A(p,iu) shown in Fig. 4, and the sohd 
curve represents an estimate of the quasiparticle dispersion using these data and Lanczos results for p 
near (0,0). The broad darkened areas represent the incoherent spectral weight in the upper and lower 
Hubbard bands. The horizontal dashed line denotes the chemical potential /^. 
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Figure 6: Particle-particle and particle-hole t-matrix equations. 
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Figure 7: (a) Singlet irreducible particle-particle vertex for zero energy transfer Fis (q, iuJm = 0) versus q 
along the (1,1) direction. Here U jt = 4 and (n) = 0.875. As the temperature decreases below J = 4:t^ /U, 
the strength of the interaction is enhanced at large momentum transfer, (b) Energy transfer dependence 
of ris(q, iwm) for q = (vr, tt), ujn = ttT, and ui„' = + Wm is shown for various temperatures. Again 
note the similarity to x(q, i'^m) in Fig. 1. 
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Figure 8: (a) Leading Q = (7r,7r), 5* = 1 particle-hole (solid circle) and the Q = 0, singlet particle- 
particle (open circle) eigenvalues versus temperature for U/t — 4 and (n) = 1. (b) Singlet d^2_^2 
eigenvalue versus temperature for U/t = 4 (solid circle) and U/t = 8 (open circle) at (n) = 0.875. 



for q ~ (tt, tt). Comparing Figs. |l] and ^ one clearly sees that the q and ujm structure of the 
interaction and x(Qi i'-^m) are similar, both reflecting the development of the antiferromagnetic 
correlations as T is reduced below J. 

Given the Monte Carlo results for the irreducible particle-particle vertex ris(p', —p',p, —p) 
in the zero energy and center-of-mass momentum channel, which we will denote by Tis{p'\p) and 
the single-particle Green's function Ga{p), the Bethe-Salpeter equation for the particle-particle 
channel is 

XMp) = -^Y.^is{p\p')G^{p')Gii-p')Mp')- (5) 
p' 

Here, as before, the sum on p' is over both p' and ujn' = (2n' + 1)ttT. Using the same Monte 
Carlo data for the four-point vertex, one can also determine the irreducible particle-hole vertex 
Ti{p\p') for a center-of-mass momentum Q = (tt, tt) and spin 1. Then, using this as a kernel, the 
Bethe-Salpeter equation for the Q = (tt, tt) particle-hole channel is 

^Mp) = E ribb')GT(y + Q)Gx(p')0Jp')- (6) 
p' 

Results showing the temperature dependence of the leading even frequency eigenvalues for the 
particle-particle (Q = 0) and particle-hole (Q = (tt, tt)) channels for a half- filled 8x8 lattice with 
U/t — 4: are plotted in Fig. ||(a). At half-filling, the dominant eigenvalue occurs in the particle- 
hole Q = (tt, tt) spin 5 = 1 channel, reflecting the development of long-range antiferromagnetic 
correlations in the half-filled ground state. The Q = (tt, tt) magnetic susceptibility varies as 

x(r)^(i-A(r))-i (7) 

and diverges as T goes to zero. The leading even frequency singlet particle-particle eigenvalue 
is associated with an eigenfunction (j){p,iuJn) which has d^^-y^ symmetry as shown in Fig. |[ 
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Figure 9: Momentum dependence of the leading singlet even-frequency particle-particle eigenfunction 
(^(p, iuJn) at uj„ — TvT for U/t — 4, (n) = 0.875 and T = 0.25t. Here, p moves on the fermi surface of the 
half-filled noninteracting system from (tt, 0) to (— 7r,0), and (?i(p, ittT) has been normalized by its value 
at p = (tTjO). The solid line denotes (cospa: —cospy)/2. Hence, near the Fermi surface 4>{p,iirT) is very 
close to the (cospx — cosp3,)/2 form. 




Figure 10: Real space structure of ris(^) for U/t = 4, (n) = 0.875 and T = 0.25t. 



Figure 11: Two-chain Hubbard lattice with inter-chain and intra-chain hopping matrix elements t± 
and t, respectively. 



Its rapid increase at low temperature reflects the fact that the doped two-hole (or two-particle) 
state is bound 

For (n) = 0.875, the Monte Carlo calculations are limited to higher temperatures because of 
the fermion determinantal sign problem. However, as seen in Fig. ^(b) for U/t = 4 and 8, the 
d^2_y2 eigenvalue increases as the temperature is lowered. This is associated with the develop- 
ment of short range antiferromagnetic correlations and the increase in Fig at large momentum 
transfer. In this context it is interesting to plot the real space Fourier transform 

(^) = ^ E ^'''■^ (q, lCUra = 0) (8) 

q 



As seen in Fig. 10, ris(£) is repulsive on site but attractive on near neighbor sites. It is this 
feature along with the large Fermi surface that leads to the d^2_y2 pairing fluctuations in the 
two-dimensional Hubbard model. 



3 Ladders 

New materials consisting of arrays of weakly coupled Cu-0 ladders have recently been synthe- 
sized [|l9| . These systems provide further opportunities to examine the physical properties of 
strongly correlated insulators and the metallic state which is reached when holes are doped into 
such insulators. Here we present a short review of the results which have been obtained for 
the two-chain Hubbard ladder using a numerical density matrix renormalization group (DMRG) 
method ||l|. 

The two-chain Cu-0 ladders in SrCuOs and La2Cu205 consist of Cu-O-Cu legs joined by 
placing an O atom between the Cu atoms of the legs, forming Cu-O-Cu rungs. In both cases, 
the ladders are also coupled to each other, but here we consider the idealized case of an isolated 
ladder. Furthermore, we approximate the ladder by the two-chain Hubbard lattice shown in 
Fig. m. Here 

(9) 

with t the hopping integral parallel to the chains, t± the hopping between the legs, and U the 
onsite Coulomb interaction. The operator cj^^, creates an electron of spin s on site («, A) with 

A = 1,2 denoting the leg and rii^xs — cl^^Ci,\s- For U = 0, H can be diagonalized in terms 
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Figure 12: Phase diagram of the U = system. 
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Figure 13: (a) Charge and (b) spin gaps, Ac and A^, versus the inverse lattice length L ^ for U/t = 8. 
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of bonding (— ) and antibonding (+) bands with Sk = 2tcoskx ±t±. A phase diagram for the 
two-chain ladder is sketched in Fig. [l^ for [/ = 0. The region above the hne t±/t = 1 — cos(7r(ri)), 
with (n) = (n^i + riii) the site fiUing, corresponds to having electrons only occupy the bonding 
band. The shaded region has both the bonding and antibonding band occupied with the bonding 
band half-occupied along the dashed curve. Thus for the half-filled case, (n) = 1, the bonding 
band is full when t±/t exceeds 2 and the non-interacting system is a simple band insulator with 
equal charge and spin gaps set by 2{t± — 2t). In the non-interacting U — Q hmit, for t^/t less 
than 2, electrons occupy both bands and the system behaves as a simple metal with vanishing 
charge and spin gaps. Here we will see that the situation is quite different in the presence of 
the Coulomb interaction U . At half-filling with [/ 7^ 0, the system remains an insulator for 
t±/t < 2, with a ratio of the spin gap to charge gap which approaches zero smoothly as t^/t 
vanishes. Furthermore, holes doped into this correlated insulator form d2.2_j,2-like pairs which 
exhibit power-law correlations over much of the shaded region, and while the charge-gap vanishes, 
a finite spin gap remains. 

The charge and spin gaps can be determined from the ground state energy Eq{N'\, Ni) of a 
configuration with iV| up-spin and down-spin electrons. The charge gap Ac is 

Ae - i [Eo{N^ + 1,N^ + 1)- Eo{N^,N^)] , (10) 

and the spin gap 

As=EoiN^ + l,N^-l)-EoiN^,N^). (11) 

Here iV| = iV| and iV| -I- iV| is fixed to give the desired filling (n) = (iVf + Ni)/ L for a 2 x L site 
ladder. Using the DMRG, calculations on ladders of varying lengths (2 x L) have been carried 
out, giving the L ^ 00 extrapolations for Ac and A^ shown in Fig. y_3| for U/t — 8 and various 
fillings. In this way, one finds that Ac and Ag are both finite for (n) = 1 but Ac vanishes for 
(n) 7^ 1. A plot of As/Ac versus t±/t for (n) = 1 is shown in Fig. |lj for several values of U/t. For 
the interacting system A^/Ac gradually approaches 1 as t±/t increases beyond 2. The difference 
(1 — As/ Ac) is one characteristic difference between the strongly correlated insulating state and 
the band insulating state. 

As noted, when holes are doped into the half-filled ladder, the charge gap vanishes. However, 
a spin gap remains which for a fixed value of U /t depends upon the ratio of i^/i and the filling. 
In the following discussion, we will set U/t — 8 and (n) = 0.875, and examine the properties of 
the system for various values of tj_/t. The dependence of the spin gap on tj_/t is illustrated in 
Fig. For t±/t > 1.7, the spin gap vanishes, consistent with the one-band Luttinger liquid 
picture suggested by Fig. |l^. For U — Q, the anti-bonding band becomes empty for (n) — 0.875 
when t±/t exceeds 1.85. As tj_/t decreases, the spin gap rises and then appears to vanish for 
t±/t ~ 0.43. This is near the point at which the bonding band is half-filled, and as discussed by 
Balents and Fisher a RNG calculation finds that the spin gap vanishes. At still lower values 
of t±/t, the RNG calculation finds a small spin gap, but we have been unable to resolve it if it 
is indeed present. 

In the spin gap regime of the doped system, the dominant two-particle correlations appear in 
the pairing channel. Figure [l^ shows the rung-rung and rung-leg pair-field correlation functions 

D^xii) ^ {Axil + l)Al{t)) (12) 
Dyx{e) = {Ay{i + £)AU^)). (13) 

Ai(0 = (cLT4i-4i4T) (14) 
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Figure 14: Ratio of the spin and charge gaps, A^/Ac, versus t±/t for U/t = A and 8 at (n) = 0.875. 
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Figure 15: Spin gap Ag versus t±/t for U /t = 8 and (n) = 0.875. 
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Figure 16: Pair field correlation functions Dxxi£) and Dyx{t) versus the lattice spacing I = \i — j\. 
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Figure 17: Power law decay of Dxx{£ = \i — j\) for different values of t±/t at U/t 
The dotted and the dashed lines decay as |i — and \i — respectively. 
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Figure 18: Comparison of the strength of the pairing correlations (open circle) with the size of the spin 
gap As (solid circle) for U/t = 8 and (n) = 0.875. The average over £ of (AiAi+f)ave has been done for 
I between 8 and 12 lattice spacings. 
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Figure 19: Magnetic structure factor S{Tr, q) versus q along the ladder for U/t = 8 and various fillings. 
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creates a singlet pair on the ith rung, while 




) 



(15) 



creates a singlet pair along a leg of the ladder at site i. The Cu-0 hopping parameters t and tj_ 
have the same sign, so that the relative phases of the pair-field correlations shown in Fig. |l^ are 
meaningful and imply that the bound pairs have _y2-\ike symmetry, as discussed by Rice et 
al. 1^ The power law decay of these pairing correlations is shown in Fig. It is also interesting 
to compare the strength of the pairing correlations with the size of the spin gap. This is done in 



These results support a picture in which the doped state of the two-leg Hubbard ladder, 
in the appropriate t±/t regime, is characterized by a finite spin gap and dx2-y2-\ike power-law 
pairing correlations. In Fig. |l9|, we show plots of the magnetic structure factor ^(Tr, q) with q 
along the ladder for various fillings. In the doped ladders, the peak in the structure factor occurs 
at an incommensurate wavevector. The exponential decay of the magnetization-magnetization 
correlation functions in the spin-gapped regime lead to Lorentzian line shapes. Thus even though 
the antiferromagnetic correlations are clearly short range, they lead to power law d,j.2__y2 pairing 
correlations. 

4 Conclusions 

The development of both the low-energy quasiparticle dispersion and the peak in the singlet 
particle-particle vertex at large momentum transfers reflects the growth of the short-range an- 
tiferromagnetic correlations as T decreases below J . Analysis of the Bethe-Salpeter equation 
shows that the leading even-frequency singlet pairing occurs in the dx2-y2 channel. Physically, 
for the large Fermi surface associated with the observed quasiparticle dispersion, a particle- 
particle vertex which increases at large momentum transfer favors dx2_y2 pairing. Note that the 
tendency for dx2_y2 pairing does not require a particularly sharp or narrow peak in Tis{p—p') for 
p — p' = (tTjTt), but rather simply weight at large momentum transfers corresponding to strong 
short-range antiferromagnetic correlations. This is clearly seen in the results for the two-leg 
ladder where the antiferromagnetic correlations exponentially decay. Thus it is the strong short- 
range antiferromagnetic correlations which lead to the formation of dx2_y2 pairing correlations 
in the Hubbard model. 

Acknowledgments 

The work on the 2D Hubbard model discussed here was carried out with N. Bulut and S.R. 
White. The work on the ladder was carried out with R. Noack and S.R. White. This work 
was supported in part by the National Science Foundation under grant No. DMR92-25027. The 
numerical computations reported in this paper were carried out at the San Diego Supercomputer 
Center. 

References 



1. J.R. Schrieffer, Solid State Commun. 92, 129 (1995) 

2. D.J. Scalapino, Phys. Reports 250, 329 (1995). 

3. D. Van Harlingen, Rev. Mod. Phys. 67, 515 (1995). 

4. D.J. Scalapino, Physica C 235-240, 107 (1994). 

5. N.E. Bickers, D.J. Scalapino and S.R. White, Phys. Rev. Lett. 62, 961 (1989). 



Fig. |T|. 



15 



6. T. Moriya, Y. Takahashi and K. Ueda, J. Phys. Soc. Japan 59, 2905 (1990). 

7. P. Monthoux and D. Pines, Phys. Rev. Lett. 69, 961 (1992). 

8. J.R. Schrieffer, J. Low Temp. Phys. 99, 397 (1995). 

9. Proceedings of the 10th Anniversary HTS Workshop on Physics, Materials and Applica- 
tions, held in Houston, March 12-16, 1996. 

10. N. Bulut, D.J. Scalapino and S.R. White, Phys. Rev. B 47, 6157 (1993). 

11. N. Bulut, D.J. Scalapino and S.R. White, Phys. Rev. B 50, 9623 (1994). 

12. R.M. Noack, S.R. White, and D.J. Scalapino, Phys. Rev. Lett. 73, 882 (1994); R.M. Noack, 
S.R. White, and D.J. Scalapino, Europhys. Lett. 30, 163 (1995). 

13. J.E. Hirsch and S. Tang, Phys. Rev. Lett. 62, 591 (1989). 

14. E. Dagotto, A. Moreo, R.L. Sugar, and D. Toussaint, Phys. Rev. B 41, 811 (1990). 

15. A. Parola, S. Sorella, M. Parrinello, and E. Tosatti, Phys. Rev. B 43, 6190 (1991). 

16. E. Dagotto, A. Moreo, F. Ortolani, J. Riera, and D.J. Scalapino, Phys. Rev. Lett. 67, 
1918 (1991); E. Dagotto, F. Ortolani, and D.J. Scalapino, Phys. Rev. B 46, 3183 (1992). 

17. N. Bulut, D.J. Scalapino, and S.R. White, Phys. Rev. B 50, 7215 (1994). 

18. A. Moreo, S. Haas, A.W. Sandvik, and E. Dagotto, Phys. Rev. B 51, 12045 (1995). 

19. Z. Hirol and M. Takano, Nature 377, 41 (1995). 

20. L. Balents and M.P.A. Fisher, preprint. 

21. T.M. Rice, M. Troyer, and H. Tsunetsugu, J. Phys. Chem. Solids 56, 1663 (1995). 



16 



